Tsallis thermostatistics for finite systems: a Hamiltonian approach 
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We show that finite systems whose Hamiltonians obey a generalized homogeneity relation rigor- 
ously follow the nonextensive thermostatistics of Tsallis. In the thermodynamical limit, however, 
our results indicate that the Boltzmann-Gibbs statistics is always recovered, regardless of the type 
of potential among interacting particles. This approach provides, moreover, a one-to-one correspon- 
dence between the generalized entropy and the Hamiltonian structure of a wide class of systems, 
revealing a possible origin for the intrinsic nonlinear features present in the Tsallis formalism that 
lead naturally to power-law behavior. Finally, we confirm these exact results through extensive 
numerical simulations of the Fermi-Pasta-Ulam chain of anharmonic oscillators. 
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I. INTRODUCTION 

Few quantities in physics play such a singular role in 
their theories as the entropy in the Boltzmann-Gibbs 
(BG) formulation of statistical mechanics: it provides, in 
an elegant and simple way, the fundamental link between 
the microscopic structure of a system and its macroscopic 
behavior. Almost two decades ago, Tsallis |Q| proposed 
the following entropy expression as a nonextensive gen- 
eralization of the Boltzmann-Gibbs (BG) formalism for 
statistical mechanics. 
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where fc is a constant, q a parameter and /(r) a probabil- 
ity distribution over the phase space variables r. Instead 
of the usual BG exponential, upon maximization with a 
fixed average energy constraint, the above entropy gives 
a power-law distribution. 
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where /3 is a temperature-like parameter and E is the 
energy of the system. In the limit q ^ 1, the above 
expressions reduce to the familiar BG forms of en- 
tropy, S = —kJ dr/(r) ln/(r), and exponential distribu- 
tion /(r) cx exp{—(3E). The consequences of this general- 
ization are manifold and far-reaching (for a recent review 
see e.g. |^), but its most notorious one is certainly the 
fact that Sq, along with other thermodynamical quanti- 
ties, are nonextensive for q ^ 1. 
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Whether such formalism can be verified experimentally 
or elucidated from any previously established theoreti- 
cal framework are obvious questions that arise naturally. 
The answer to the first one seems to be affirmative, as 
there are currently numerous references to systems that 
are better described within the generalized approach of 
Tsallis than with the traditional BG formalism |||, ^, the 
underlying argument being usually based on a best fit 
to either numerical or experimental data by choosing an 
appropriate value for the parameter q. This is, by far, 
the most frequent approach towards Tsallis' distribution 
and is surely not conclusive. Its justification from first 
principles, although already addressed at different levels 
in previous studies |Q, ||, |^ , is nevertheless controversial 
and still represents an open research issue. For example, 
it has been recently suggested that systems in contact 
with finite heat baths should follow the thermostatistics 
of Tsallis 1^. This derivation, however, relies on a partic- 
ular ansatz for the density of states of the bath and lacks 
a more direct connection with the entropy, as pointed 
out in Ref. 0. Also, it has been shown in Ref. |g] that 
the Tsallis statistics, together with the biased averaging 
scheme, can be mapped into the conventional Boltzmann- 
Gibbs statistics by a redefinition of variables that results 
from the scaling properties of the Tsallis entropy. 

In a recent study |9j, a derivation of the generalized 
canonical distribution is presented from first principles 
statistical mechanics. It is shown that the particular fea- 
tures of a macroscopic subunit of the canonical system, 
namely, the heat bath, determines the nonextensive sig- 
nature of its thermostatistics. More precisely, it is exactly 
demonstrated that if one specifies the heat bath to satisfy 
the relation 
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where /3 oc 1/ (qkT) is a temperature-like parameter and 
T is the thermodynamic temperature, the form of the 



2 



distribution Eq. (Q) is recovered Equation ^ is es- 
sentially equivalent to Eq. (||). However, it reveals a di- 
rect connection between the finite aspect of the many- 
particle system and the generalized g-statistics It 
is analogous to state that, if the condition of an infinite 
heat bath capacity is violated, the resulting canonical dis- 
tribution can no longer be of the exponential type and 
therefore should not follow the traditional BG thermo- 
statistics. Inspired by these results, we propose here a 
theoretical approach for the thermostatistics of Tsallis 
that is entirely based on standard methods of statistical 
mechanics. Subsequently, we will not only recover the 
previous observation that an adequate physical setting 
for the Tsallis formaHsm should be found in the physics 
of finite systems, but also derive a novel and exact cor- 
respondence between the Hamiltonian structure of a sys- 
tem and its closed-form g-distribution, supporting our 
findings through a specific numerical experiment. 



II. THEORETICAL FRAMEWORK 

We start by considering, in a shell of constant energy, 
a system whose Hamiltonian can be written as a sum of 
two parts, viz. 



and computing 
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where r 
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(ri, . . . ,r2Ar), with 
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5 ' ni+n2 



+„2-i-i,--- j^ni+n2+n3) and so on. The fact that 
Tsallis distribution is a power-law instead of exponen- 
tial strongly suggests us to look for scale-invariant forms 
of Hamiltonians Furthermore, since scale- invariant 
Hamiltonians constitute a particular case of homoge- 
neous functions our approach here is to show that, 
if H2 satisfies a generalized homogeneity relation of the 
type 



Ai/2(r2, rg, . . . ) = iJ2(Ai/'^^r2, A^/'^^rg, 
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where X,a2,a3, ■ ■ ■ are non-null real constants, then the 
correct statistics for i/i(ri) is the one proposed by Tsal- 
lis. The foregoing derivation is based on a simple scaling 
argument, but we shall draw parallels to Ref. j9| when- 
ever appropriate. 

The structure function (density of states) for H2 at the 
energy level E is given by 



n2iE)^ J S{H2{y2,V^, . . .) - E)AV2. 
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where dV2 = d"^r2d"^r3 ... is the volume element in the 
subspace spanned by (r2, rg, . . . ). For systems satisfying 
equation (||) , this function can be evaluated taking A > 



f^2(A-Bo) = J S{H2iv2,r3,...)-XEo)dV2 

= jJs{X-'H2{r2,r3,...)-EQ)dV2 

= J f S{H2iX'^r2,X'^r3,...)-Eo)dV2 
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= Ai/(^-i)l]2(i?o), 
where we define 

1 = ^ 
a - 1 ~ . ^ a. 
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1=2,3, 



and utilize the notation r2 = A ^/'^^r2, rg — A ^^"■^r^, 
etc., and dV^ = d""r^d""r(, . . . Hence, if O2 is defined 
at a value Eq, it is also defined at every E = XEq, with 
A > 0. We can then write 



(9) 



and express the canonical distribution law over the phase 
space of Hi as 



fin) 
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where H is the total energy of the joint system composed 
by Hi and H2, and Q{H) is its structure function, 

n{H)^ J 5{Hi{vi) + H2{y2,V3, . . .) ~ H)dVidV2, 

(11) 

where dVi and dV2 are the infinitesimal volume elements 
of the phase spaces of Hi and H2 , respectively. Compar- 
ing Eq. ( [lo| ) with the distribution in the form of Eq. (||), 
we get the following relation between 5, /3 and H: 
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Notice that one could reach exactly the same result us- 
ing the methodology proposed in Ref. |^, i.e. by evalu- 
ating I3{E) = n'2{E)/n2{E) &t E ^ H, calculating q - 1 
through Eq. (^) and inserting these quantities back in 
Eq. d). 

As already mentioned, previous studies have shown 
that the distribution of Tsallis Eq. (||) is compatible 
with some anomalous "canonical" configurations where 
the heat bath is finite [Q| or composes a peculiar type of 
extended phase-space dynamics ||l^. In our approach. 
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the observation of TsalUs distribution simply reflects the 
finite size of a thermal environment with the property 
the thermodynamical Hmit corresponding to g 1 
in Eq. (^. We emphasize that, although similar conclu- 
sions could be drawn from Refs. |^, the theoretical 
framework introduced here permits us to put forward a 
rigorous realization of the g-thermostatistics: it stems 
from the weak coupling of a system to a "heat bath" 
whose Hamiltonian is a homogeneous function of its co- 
ordinates, the value of q being completely determined by 
its degree of homogeneity, Eq. This provides also 

a direct correspondence between the parameter q and 
the Hamiltonian structure through geometrical elements 
of its phase space, viz. the surfaces of constant energy 
H2 = E- 

As a specific application of the above results, we in- 
vestigate the form of the momenta distribution law for a 
classical iV-body problem in d-dimensions. The Hamil- 
tonian of such a system can be written as 



ff(p,q) 



1 ^ 

i=l 
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where we define Hi = ^ J2tLi Pi > P» 



[Pii 



,qw),(i3) 

■ • ,Pid) is 

the linear momentum vector of an arbitrary particle i 
(hence the number of degrees of freedom of the system 1 
is ni = Nd), and H2 (the "bath") is due to a homogeneous 
potential V{q) of degree a, i.e., XV{c[) = F(Ai/"q) with 
q = (qi, • ■ • , (In)- At this point, we emphasize that the 
distinction between "system" and "bath" is merely formal 
and does not necessarily involve a physical boundary. It 
relies solely on the fact that we can decompose the total 
Hamiltonian in two parts By making the correspon- 
dences ri = (pi, . . . ,pn), r2 = (qi, . . . ,qA'), ni = Nd, 
712 = Nd, ai — 2 and 02 — a, the homogeneity relation 
(^ is satisfied. From Eq. (10) it then follows that 
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where the nonextensivity measure q is given by 
1 Nd Nd 



Nd 

= 1^0 = 

la Nd-a 
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It is often argued that the range of the forces should play 
a fundamental role in deciding between the BG or Tsallis 
formaHsms to describe the thermostatistics of an iV-body 
system ||^. For example, the scaling properties of the 
one-dimensional Ising model with long-range interaction 
has been investigated analytically IQ and numerically 
1^ in the context of TsaUis thermostatistics, whereas 
in Ref. ||l^ a rigorous approach was adopted to study 
the nonextensivity of a more general class of long-range 
systems in the thermodynamic limit (see below). Recall 
that, for a d-dimensional system, an interaction is said 
to be long-ranged if —d < a < 0. Within this regime. 



the thermostatistics of Tsallis is expected to apply, while 
for a < —d the system should follow the standard BG 
behavior ||2j. This conjecture is not confirmed by the re- 
sults of the problem at hand. Indeed, Eq. (|l4) is consis- 
tent with the generalized g-distribution Eq. ^ no matter 
what the value of a is, as long as it is non-null and N 
is finite. In the Hmit TV ^ 00, however, we always get 
q ^ 1, with the value of a determining the shape of the 
curve q = q{N). If a > 0, g approaches the value 1 from 
above, while for a < the value of q is always less than 
1. Therefore, for (ergodic) classical systems with N par- 
ticles interacting through a homogeneous potential, the 
equilibrium distribution of momenta always goes to the 
Boltzmann distribution, /(ri) oc exp[— /?i/i(ri)], when 
TV — > 00. This observation should be confronted with the 
recent results of Vollmayr-Lee and Luijten ||l^, who in- 
vestigated the nonextensivity of long-range (therein "non- 
integrable") systems with algebraically decaying interac- 
tions through a rigorous Kac-potential technique. Con- 
trary to the trend established by the practitioners of Tsal- 
lis' formaHsm, those authors argue that it is possible to 
obtain the nonextensive scaling relations of TsaUis with- 
out resorting to an a priori g-statistics, the Boltzmann- 
Gibbs prescription (g = 1) being sufficient for describing 
long-range systems of the type above. Even though our 
findings embody partially the same message (we are not 
yet concerned about scaling relations), there are some 
caveats that prevent their results from being straightfor- 
wardly applicable to our problem: neither a system-size 
regulator for the energy nor a cutoff function is present in 
our treatment. This is immediately in contrast with their 
observation that the "bulk" thermodynamics strongly de- 
pends on the functional form of the regulator. Moreover, 
by not addressing the distribution function expHcitly at 
finite system sizes, that work has very little in common 
with the most interesting part of our study, which might 
in fact explain some observations of the g-distribution. 
Notwithstanding these differences, we beheve that an in- 
vestigation of the scaling properties of the system stud- 
ied here would elucidate from a different perspective the 
connection of TsaUis thermostatistics with nonextensiv- 
ity and is certainly a very welcome endeavor. 

It is important to stress here that the essential fea- 
ture determining the canonical distribution is the geom- 
etry of the phase space region that is effectively visited 
by the system. In a previous work by Latora et. al. 
Jl8| , the dynamics of a classical system of N spins with 
infinitely long-range interaction is investigated through 
numerical simulations, and the results indicate that if 
the thermodynamic limit {N 00) is taken before the 
infinite-time limit {t 00), the system does not relax 
to the Boltzmann-Gibbs equUibrium. Instead, it dis- 
plays anomalous behavior characterized by stable non- 
Gaussian velocity distributions and dynamical correla- 
tion in phase space. This might be due to the appear- 
ance of metastable state regions that have a fractal na- 
ture with low dimension. In our theoretical approach, 
however, it is assumed that the infinite-time limit is 



taken before the thermodynamic limit. As a consequence, 
metastable or quasi-stationary states like the ones ob- 
served by Latora et al. [IS) with a particular long- 
range Hamiltonian system cannot be predicted within 
the framework of our methodology. Whether this type 
of dynamical behavior can be generally and adequately 
described in term of the nonextensive thermostatistics of 
Tsallis still represents an open question of great scientific 
interest. 



III. NUMERICAL EXPERIMENTS 

In order to corroborate our method, we investigate 
through numerical simulation the statistical properties 
of a linear chain of anharmonic oscillators. Besides the 
kinetic term, the Hamiltonian includes both on-site and 
nearest-neighbors quartic potentials, i.e. 



N 



N 



N 



H 



{Qi+i - Qt) 



(16) 



The choice of this system is inspired by the so-called 
Fermi-Pasta-Ulam (FPU) problem, originally devised to 
test whether statistical mechanics is capable or not to 
describe dynamical systems with a small number of par- 
ticles ||l^. From Eq. (16), we obtain the equations of 



motion and integrate them numerically together with the 
following set of initial conditions: 



fc=l 



2Trik 
cos I -I- ^pk 



,N (17) 



where is a random number within [0, 27r). Undoubt- 
edly, a rigorous analysis concerning the ergodicity of this 
dynamical system would be advisable before adopting the 
FPU chain as a plausible case study. This represents a 
formidable task, even for such a simple problem pO[ . For 
our practical purposes, it suffices, however, to test if the 
system displays equipartition among its linear momen- 
tum degrees of freedom, since this is one of the main sig- 
natures of ergodic systems. Indeed, one can show from 
the so-called Birkhoff-Khinchin ergodic theorem that, for 
(almost) all trajectories |^ , 

/ dH\ , , ViE) , , 

where r(i?) = /^(^j^^ dV is the volume of the phase 

space with B <E, {A)(t) = (1/t) /q AtA{t) denotes the 
the usual time average of an observable A, and 9 stands 
for the absolute temperature of the whole system, Q = 
4E/{3N) (cf. 11). We then follow the time evolution of 
the quantities 



1 

e 



dH 

' dpi 



(t), t = 1, 



(19) 




log (H-H^) 



FIG. 1: Logarithmic plot of the rescaled distribution 
Hl-^'^f{H - Hi) as a function of the transformed variable 
H — Hi for N = 8 (circles), 16 (triangles up), 32 (squares), 
and 64 (triangles down) anharmonic oscillators. The solid 
straight lines are the best fit to the simulation data of the 
expected power-law behavior Eq. ([T^). The slopes are 1.0068 
(1.0), 3.07 (3.0), 7.21 (7.0), and 15.27 (15.0) for = 8, 16, 32, 
and 64, respectively, and the numbers in parentheses indicate 
the expected values obtained from Eq. (|l5[). The departure 
from the power-law behavior at the extremes of the curves is 
due to finite-time sampling. 



to check if they approach unity as t increases. From the 
results of our simulations with different values of N and 
several sets of initial conditions, we observe in all cases 
the asymptotic behavior, Qi{t) — > 1 as t ^ oo. This pro- 
cedure also indicates a good estimate for the relaxation 
time of the system, r « 10^, so we shall consider our sta- 
tistical data only for t > 10^, with a typical observation 
time in the range 10 — 10*, after thermalization. 

In Fig. 1 we show the logarithmic plot of the dis- 
tribution /(pi,P2, • ■ ■ ) against the transformed variable 
H ~ Hi for systems with A^ = 8, 16, 32, and 64 oscil- 
lators. Based on the above result, we assume ergodicity 
and compute the distribution of momenta from the fiuc- 
tuations in time of Hi through the relation 



/(iIl)«i/}^/')-V(Pl,P2,...), 



(20) 



where the h[^^^''~^ factor accounts for the degeneracy 
of the momenta consistent with the magnitude of Hi (cf. 
lUl). Indeed, we observe in all cases that the fiuctuations 
in Hi follow very closely the prescribed power-law behav- 
ior Eq. (|lj) , with exponents given by Eq. ( |l^ ) . These re- 
sults, therefore, provide clear evidence for the validity of 
our dynamical approach to the generalized thermostatis- 
tics. 



IV. CONCLUSION 

In conclusion, we have shown that the generalized 
formalism of Tsallis can be applied to homogeneous 
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Hamiltonian systems to engender an adequate theoret- 
ical framework for the statistical mechanics of finite sys- 
tems. Of course, we do not expect that our approach can 
explain the whole spectrum of problems in which Tsal- 
lis statistics can be applied. However, our exact results 
clearly indicate that, as far as homogeneous Hamiltonian 
systems are concerned, the range of the interacting po- 
tential should play no role in the equilibrium statistical 
properties of a system in the thermodynamic limit | p3[ . 
Under these conditions, the conventional BG thermo- 
statistics remains valid and general, i.e., for the specific 
class of homogeneous Hamiltonians investigated here, the 



thermodynamic limit {N co) leads always to BG dis- 
tributions. 
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